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The 18 Ne(a,p) 21 Na reaction is one key for the break-out from the hot CNO-cycles to the rp- 
process. Recent papers have provided reaction rate factors Na < crv > which are discrepant by at 
least one order of magnitude. The compatibility of the latest experimental results is tested, and a 
partial explanation for the discrepant Na < av > is given. A new rate factor is derived from the 
combined analysis of all available data. The new rate factor is located slightly below the higher rate 
factor by Matic et al. at low temperatures and significantly below at higher temperatures whereas 
it is about a factor of five higher than the lower rate factor recently published by Salter et al. 

PACS numbers: 25.60.-t,25.55.-e,26.30.-k 



I. INTRODUCTION 

The reaction rate of the 18 Ne(a,p) 21 Na reaction pro- 
vides a route from hot CNO-cycles to the NeNa and MgAl 
cycles and finally to the rp-process at typical tempera- 
tures of e.g. about 1 — 2 GK (Tg = 1 — 2) in X-ray bursters 
(XRB) pj. It is expected that this reaction is the domi- 
nating route in the low temperature range Q. An alter- 
native route from hot CNO-cycles to the rp-process may 
be the 15 0(a,7) 19 Ne reaction. 

The relatively high temperatures correspond to most 
effective energies of about 1.3 to 2.1 McV for the 
18 Ne(ct,p) 21 Na reaction which are experimentally well ac- 
cessible. However, experiments remain very difficult be- 
cause of the short-living 18 Ne nucleus (X\/ 2 = 1.67 s) 
and the limited intensity of radioactive beam facilities. 
Thus, besides the direct approach of measuring the 
18 Ne(a,p) 21 Na reaction cross section 0, the reverse 
21 Na(p,ct) 18 Ne reaction has been studied very recently 
[f| and in an earlier unpublished experiment @ , and the 
resonance energies have been determined from various 
transfer experiments populating states in the compound 
22 Mg nucleus [Mil- 

The focus of the present paper is the comparison of 
the latest experiments by Groombridge et al. (hereafter: 
GRO) Q, Salter et al. (SAL) @, Chae et al. (CHA) [T^], 
and Matic et al. (MAT) [ll|- The earlier direct data of 
[H have been improved and extended by the same group 
leading to the GRO data. The SAL data are the only 
published data for the inverse reaction; a brief compar- 
ison to the unpublished data measured at Argonne Na- 
tional Laboratory (ANL) is also provided. The MAT 
transfer data have by far the best energy resolution which 
is essential for a precise determination of the resonance 



energies. Additional measurements of angular distribu- 
tions in CHA lead to a new spin assignment only in few 
cases (see Table III of CHA). 

The reaction rate factor N A < crv > for the 
18 Ne(a,p) 21 Na reaction is given by the sum over the con- 
tributing resonances: 



cm 3 s _1 mol~ 



1.54 x 10 11 GuT 9 )- 3/2 

^(^7). x exp(-11.605^/T 9 )(l) 



with the reduced mass [i in units of amu, the resonance 
energies Ei in MeV, and the resonance strengths (0^7); in 
MeV. In general, resonance energies are given as E in the 
center-of-mass (cm.) system without index; excitation 
energies are given as E* in this paper. 

The resonance strength (017) for the 18 Ne(a,p) 21 Na re- 
action is given by 



cj 7 = (2J+1) 



r r 
r 



(2) 
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with the resonance spin J, the partial widths T a and F p , 
and the total width L = T a + T p + P 7 w T a + T p . In 
most cases it can be expected that T a <C T p . and thus 
(td"f) ~ (2 J+l) T a . The application of the simple formula 
for narrow resonances in Eq. (JXJ) is justified because the 
resonance widths L are much smaller than the resonance 
energies E [l2l ]. 

In the following we first briefly review the various ex- 
perimental approaches and discuss the resulting uncer- 
tainties in the determination of the reaction rate factor 
Na < <jv >. Next we check whether the experimental 
results of GRO, MAT, CHA, SAL, and ANL are com- 
patible with each other. Finally, the reaction rate fac- 
tors Na < crv > of the different studies are compared. 
Note that Na < crv > of different studies may differ not 
only from discrepant resonance energies and resonance 
strengths or cross sections, but also from a different num- 
ber of considered resonances in Eq. (JXJ) . 
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II. BRIEF REVIEW OF AVAILABLE DATA 
A. Transfer data 

Various transfer experiments have been performed in 
the last decade to study properties of the compound nu- 
cleus 22 Mg A detailed comparison of the results 
is provided in MAT and is not repeated. Here we briefly 
summarize the MAT results and some modifications re- 
sulting from the CHA data. 

Transfer data are able to provide excitation energies 
E* and spin and parity J 11 of states in 22 Mg. However, 
from the transfer data it is not possible to determine res- 
onance strengths W7 which are the second ingredient for 
the calculation of the rate factor Na < crv > in Eq. (JlJ . 

1. Matic et al. 

In the MAT approach the 24 Mg(p,t) 22 Mg reaction is 
used to populate excited states in the compound nu- 
cleus 22 Mg at proton energies of slightly below 100 MeV. 
The experiment has been performed using the Grand 
Raiden spectrometer at RCNP, Osaka. The excellent 
energy resolution of about 13keV allows a precise de- 
termination of excitation energies E* which enter expo- 
nentially into Na < crv > in Eq. flTJ via E* = E + S a 
(with the separation energy of the a particle in 22 Mg of 
S a = 8.142 MeV) and are thus the main source of uncer- 
tainties. (S a = 8.142 MeV is taken from the new Audi 
and Meng compilation [l3j]; the small difference to the 
earlier result of S a = 8.139 MeV [lH does practically not 
affect the rate factor Na < <tv > in the relevant temper- 
ature range around Tg = 1 — 2.) 

In addition to the excitation energies E* , the total 
widths r can be determined from these data by fitting the 
observed peak widths r o bs- Most of the observed states 
are much broader than the experimental resolution, and 
thus the required unfolding procedure leads only to mi- 
nor additional uncertainties for the derived width V. The 
results arc listed in Table [TJ As can be seen from Table Q] 
practically all resonances fulfill the criterion of T/E < 0.1 
which is often used as definition for narrow resonances 
(although also more stringent definitions for narrow res- 
onances can be found in literature). As we will show in 
Sect. Ill A 21 the simple formula for narrow resonances in 
Eq. (QJ provides the reaction rate factor Na < <jv > for 
the 18 Ne(a,p) 21 Na reaction with sufficient accuracy. In 
this sense the resonances in Table U can be considered 
generally as narrow resonances. 

The resonance strengths ui^f in MAT have gener- 
ally been calculated by the following procedure. In 
a first step spin and parity J 7 * for the states seen in 
the 24 Mg(p,i) 22 Mg experiment were tentatively assigned 
from the mirror nucleus 22 Ne. Note that the level scheme 
for the stable mirror nucleus 22 Ne is well-established up 
to relatively high excitation energies. In some cases also 
theoretical predictions from the shell model have been 



TABLE I: Excitation energy E* , resonance energy E, spin 
and parity ,E , total width V, and resonance strength 0)7 
for excited states in 22 Mg from the 24 Mg(p,t) 22 Mg experi- 
ment in |llf | . Later revisions for individual states are marked 
by "* GRO " and "* CHA " ; these revisions are based on the re- 
placement of the experimental resonance strengths of GRO 
by calculated resonance strengths and on revised (but still 
tentative) spin assig nments by CHA (see also Table |H}. 
The finally recommended strengths will be slightly lower by 
a factor of 0.55 (see discussion in Sect. IIII "Bl and ITV)| . 
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used [15j . Next, upper limits for T a were calculated from 
the single-particle Wigner limit using the tentative spin 
assignments (see Table VII in MAT, partly repeated in 
Table H| in this work). Finally, the Wi gner limit has been 
scaled by carefully chosen reduced widths. Whenever 
possible, the reduced a-widths were taken from experi- 
mental data obtained for the stable mirror nucleus 22 Ne 
by a-transfer on the mirror target 18 O. In the remain- 
ing cases where no experimental information is available, 
the reduced a-widths were estimated by simple but rea- 
sonable theoretical assumptions. A detailed discussion is 
given in Sect. V in MAT. 

In the few cases where experimental data are available 
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from the GRO study, experimental resonance strengths 
were used by MAT in their calculation of Na < crv >. 
Levels from threshold up to an excitation energy of about 
13MeV (corresponding to E ps 5MeV) are taken into 
account in MAT (see their Table VII), and thus the re- 
sulting rate factor Na < uv > is well determined over 
a broad temperature range including the astrophysically 
most relevant range of Tg = 1 — 2. 

For an independent comparison of the various studies 
of the 18 Ne(a,p) 21 Na reaction, we replace the experimen- 
tal resonance strengths from the GRO data by calcula- 
tions similar to the other resonance strengths in MAT 
(see upper part of Table fll] these data will be referenced 
as MAT-th: "th" for theoretical strengths only). A fur- 
ther explanation for this replacement will become visible 
later in the comparison of the MAT data to the GRO 
data and SAL data (see Sect. IIII Cl) . 



TABLE II: Revisions for resonance strengths W7, compared 
to Table VII in MAT and Table [I] The upper five lines are re- 
calculations of ix>7 instead of the resonance strengths adopted 
from GRO (here ij7mat is identical to <^>7gro)- The lower five 
lines result from new spin assignments in CHA. For compar- 
ison, the original strengths lj^/mat from MAT are also listed. 



gies between 2 and 3 MeV because of the high resonance 
strengths taken from the GRO data. 
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Under the realistic assumption that the total width L 
is dominated by the proton partial width L p (and thus 
CJ7 « u)T a and L pa r p ), the reaction cross section a(E) 
as a function of energy can be calculated as a sum over 
Breit-Wigner resonances: a(E) = J2i a BW,i(E) with 
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FIG. 1: (Color online) Reaction cross section of the 
18 Ne(a,p) 21 Na reaction, calculated from total widths F (see 
Table HJ and from resonance strengths ury. The full black line 
uses the MAT strengths with all revisions in Table [IT] (refer- 
ence cross section <T re f.). The green short-dashed line refers to 
the original MAT strengths (from their Table VII). The long- 
dashed magenta line shows the five GRO resonances only. 
The black dotted line (MAT-th) uses the MAT strengths but 
replaces the experimental GRO strengths by the calculated 
strengths (see first five lines in Table [TTJ) . MAT-th is almost 
identical to MAT (green short-dashed) at low energies and 
to oVef. (full black) at higher energies; thus, the black dot- 
ted MAT-th curve becomes only visible at energies around 
E w 1.8 MeV. The experimental data points have been deter- 
mined by the inverse 21 Na(p,a) 18 Ne reaction and represent 
the ground state contribution only. Further discussion see 
text. 
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The result is shown in Fig. [T] (MAT-th, thin dotted black 
line). The upper limit of the total width T has been 
used for the three resonances at E* = 9.157, 9.482, and 
9.542 MeV. The resulting uncertainty from these upper 
limits for the reaction rate factor Na < crv > remains 
negligible because the resonances are narrow in any case. 

This calculation of the cross section a(E) enables a 
more detailed comparison of the data from transfer and 
from the study of the inverse reaction (see Sect. IIII BP . 
The original strengths of MAT (using the experimental 
GRO strengths where available) lead to the green short- 
dashed curve in Fig. Q] which is much higher at ener- 



2. Chae et al. 

The CHA data are also based on a study of the 
24 Mg(p,f) 22 Mg reaction. The experiment has been per- 
formed at lower energies slightly above 40 MeV at the 
Holifield Radioactive Ion Beam Facility at Oak Ridge. 
The achieved energy resolution does not reach the MAT 
data, but the larger angular range under study allows an 
improved but still tentative assignment of spin and parity 
J 77 of several states in 22 Mg. 

We have taken the tentative assignments from the 
CHA data and have recalculated the resonance strengths 
W7 for these states using the adopted energies of the high- 
resolution MAT data. The results are listed in Table [TT1 
(lower part). The new spin assignments lead to smaller 
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spins J for some low-lying resonances and consequently 
to larger resonance strengths wy because of the reduced 
centrifugal barrier. However, the rate factor Na < av > 
is given as the sum over many resonances, and the influ- 
ence of the new spin assignments on the total rate factor 
remains below a factor of two over the astrophysically 
relevant temperature range of Tg = 1 — 2. 

The calculated resonance strengths us-f from Table VII 
in MAT with all modifications from Table [TT] are taken 
as a reference in this work to calculate <7 re f . {E) and the 
rate factor Na < av > ro f.. For simplicity, the rate factor 
Na < av > has been calculated from the approximation 
in Eq. ((T|). Thus, it will also be possible to modify our 
rate factor without much efforts as soon as improved ex- 
perimental data for any of the 18 Nc(a,p) 21 Na resonances 
will become available. The result is shown as the full 
black line in Fig. [TJ Note that the final recommendation 
for a(E) and Na < av > from the combined analysis of 
all available experimental data will be almost a factor of 
2 lower (see Sect. HV)) . 

We have tested the validity of the simple rate formula 
for narrow resonances in Eq. (flj by a numerical integra- 
tion of the cr ro f .(£') curve, and it is found that the devi- 
ation between the numerical integration and the simple 
formula in Eq. (TTJ) remains below 10 % over the full tem- 
perature range under study (Tg = 0.25 — 3) and far below 
5 % for the most relevant range of Tg = 1 — 2. 

It is obvious that the calculated cross section a{E) 
and reaction rate factor Na < av > depend on the spin 
and parity of the considered resonances. A detailed 
study of the corresponding uncertainties has already been 
given in MAT using random spin assignments, and it was 
concluded that Na < av > does not change by more 
than one order of magnitude. Here we provide reduction 
factors r s p .(L = 0)/r s . p . (L) for angular momenta < 
L < 4 in the energy range between 1 and 4MeV (see 
Table IIII[) where the a single-particle widths T s , p , have 
been calculated similar to MAT. The largest reduction 
factors are found for large angular momenta L at low 
energies where the influence of the centrifugal barrier is 
most important. Here an increase of L by 1 leads to a 
reduction of a factor of 10 whereas at higher energies and 
small L the increase of 1 leads to a reduction of less than 
a factor of 2. 

TABLE III: Reduction factors T S _ P .(L = 0)/T s . p .(L) for an- 
gular momenta < L < 4 in the energy range between 1 and 
4MeV. 
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The assignment of spin and parity J 7 ' of the states in 
22 Mg (e.g. in MAT) is tentative and has mainly been 
made from the mirror nucleus 22 Ne. Even if some of the 
individual assignments may be incorrect, the distribution 



of spins should be correct because the J ir assignments in 
22 Ne are mainly firm and not tentative. The cross sec- 
tion a(E) of the 18 Ne(a,p) 21 Na reaction is composed of 
overlapping resonances, and the calculation of the rate 
factor Na < av > requires folding with the Maxwell- 
Boltzmann distribution. Thus, at any astrophysically 
relevant temperature Tg > 1 the rate factor is defined 
by a sum over several contributing resonances, and the 
influence of a modified spin assignment of an individual 
resonance remains very limited, in particular, if the spin 
distribution is kept; i.e., the modification of J* of one 
resonance (e.g. increase of L and reduced r a and wy) is 
compensated by a similar modification of J x of another 
resonance (decrease of L and enhanced T a and ujj). 

As we will see later, the modified spin assignments of 
the CHA experiment (mainly smaller J n than used by 
MAT, sec Tablc[TTl lower part) will lead to a marginal en- 
hancement of the rate factor Na < av > around Tg « 1 
by less than a factor of 2. In combination with the above 
arguments, it seems thus reasonable that the uncertainty 
of the spin assignments for the calculated rate factor 
Na < av > does not exceed a factor of two in the full 
temperature range under study. 



B. Direct data 

The GRO experiment has been performed at the Ra- 
dioactive Ion Beam facility at Louvain-la-Neuve. They 
have measured the 18 Ne(a,p) 21 Mg reaction directly us- 
ing an extended 4 He gas target and a 18 Ne beam. The 
chosen detection technique allowed the reconstruction of 
the interaction vertex in the extended gas target, and 
together with the measured proton energy it was pos- 
sible to determine the energies and resonance strengths 
of 8 resonances and their main decay branch into the pi 
channel (po corresponds to the 21 Na ground state, p\ to 
the first excited state, and so on). The data cover the 
energy range from about 1.7 to 2.9 MeV. The rate factor 
Na < av > was calculated from these 8 resonances; this 
Na < av > is considered as a lower limit by GRO be- 
cause of the missing resonances outside the studied en- 
ergy region and because of perhaps missed weak reso- 
nances or missed weak branches of observed resonances 
inside the studied energy region. 

The cross section a(E) as a function of energy is calcu- 
lated from 5 adopted resonances (sec discussion below in 
Sect. IIII A)) using excitation energies E* from MAT, total 
widths r from Table [H and resonance strengths ojj from 
GRO. The result is shown in Fig. [T] with a long-dashed 
magenta line. Obviously, in the energy region between 2 
and 3 MeV it is close to the MAT result because MAT 
have used resonance strengths from GRO, but it is much 
higher than the reference calculation of this work. 

Finally, it is interesting to note that the summed po 
resonance strength to the ground state is about 42 % of 
the total summed strength in the GRO data. This num- 
ber will be relevant for comparison with the SAL and 
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ANL data for the inverse 21 Na(p,a) 18 Ne reaction. 



C. Data for the inverse reaction 

The latest experiment by SAL and the unpublished 
ANL experiment have used the reverse 21 Na(p,a) 18 Ne 
reaction in inverse kinematics with a radioactive 21 Na 
beam and a solid CH2 target. An average cross section 
a at the energy E c r is determined from the measured a 
yield: 



HE eff ) 



1 



2AE 



lE e ff-AE 



a(E) dE 



(4) 



where E e g is the energy in the center of the CH2 target, 
and the total energy loss in the target is given by 2 AE. 
In the energy range under study the 21 Na(p,a) 18 Ne cross 
section populates mainly the 18 Ne ground state (no event 
to the first excited state in 18 Ne is observed by SAL). 
Thus, the measured 21 Na(p,a) 18 Ne cross section can be 
converted to the 18 Ne(a,po) 21 Na g . s . cross section using 
detailed balance. 

Because of the relatively thick targets that are used 
in the reverse reaction experiments, the energy of each 
data point is not very well defined. This leads to major 
uncertainties in the calculation of the reaction rate factor 
Na < <jv >. For example, for the lowest data point of 
SAL one finds a variation of the astrophysical S-factor 
of a factor of 5 if the given cross section is converted to 
the S-factor at the upper or lower energy limit of the 
data point. Therefore a different way for the comparison 
of rate factors from the various experimental techniques 
has been chosen in this work (see Fig. [5] and discussion 
of Eq. © in Sect. HUB}. 



1. Salter et al. 

The SAL experiment used the ISAC II facility at TRI- 
UMF. The energy range under study by SAL for the (a,p) 
reaction (from about 1 McV up to about 2.6 McV) ex- 
tends the energy range of GRO down to lower energies 
and reaches the Gamow window of the 18 Ne(a,p) 21 Na 
reaction for the first time. The SAL data are shown in 
Fig. [T]as blue points. 

From the average cross sections of the (a,po) reaction 
the rate factor Na < crv > is calculated in SAL using 
the exp2rate code [l6|]. Because of the missing contri- 
butions of the (a,pi^o) channels this Na < av > is also 
considered as a lower limit in SAL. These missing con- 
tributions are estimated in SAL from Hauser-Feshbach 
(HF) calculations by T. Rauscher. The calculations in- 
dicate that the total rate factor Na < ov > is about 
a factor of three larger than the measured ground-state 
contribution. 

The validity of HF calculations is somewhat uncertain 
in the present case because of the relatively small level 



density in the compound nucleus 22 Mg. However, this 
uncertainty mainly influences the absolute value of the 
calculated cross section which depends on the number of 
states (or resonances) at the energy under study. The cal- 
culation of the decay branches in the proton channel to to 
the ground state and to excited states in 21 Na should be 
less affected because it is dominated by the transmission 
coefficients. It can also be checked "by hand" by calcu- 
lating single-particle limits for the proton decay of states 
in 22 Mg. It is found that e.g. the decay of a hypothetical 
0+ state in 22 Mg at E* = 9.942 MeV (E = 1.8 MeV, i.e. 
in the center of the analyzed energy region of Fig. [T]) by 
proton emission proceeds by 31 % to the ground state of 
21 Na and by 69 % to excited states. This confirms the HF 
approach for the decay branch. A further confirmation 
is obtained from comparison with the GRO data. Here 
it is found experimentally that the ground state branch 
contributes with 42 % to the summed strength. 



2. ANL data 

A similar experiment has also been performed at 
the Argonne National Laboratory (ANL). Unfortunately, 
these data have never been published and can be found 
in the ANL Annual Report only [f| . The data (extracted 
from the figure in are shown in Fig. Q] as red squares. 

The ANL data point at the highest energy agrees well 
with SAL. A second data point at lower energies is about 
a factor of two higher, and the error bars of ANL and SAL 
are close to overlap. Three upper limits have been deter- 
mined at lower energies; these upper limits are slightly 
higher than the SAL data points. In total, the unpub- 
lished ANL data are in reasonable agreement with the 
SAL data. 

As usual, it is very difficult to estimate the reliability 
of unpublished data, and thus a publication of the ANL 
data would be very helpful. The ANL data do not enter 
directly into the recommended rate (see Sect. II Vp. but the 
reasonable agreement with SAL strengthens the validity 
of the experimental data by SAL. 



III. COMPATIBILITY OF RESULTS FROM 
VARIOUS EXPERIMENTAL TECHNIQUES 

In this section we will first analyze whether the ex- 
perimental data of the various experimental approaches 
are compatible with each other. In a second step the 
uncertainties of additional ingredients for the calculation 
of the rate factor Na < crv > will be studied. These 
are in particular calculated resonance strengths W7 for 
the transfer experiments and the theoretically estimated 
ground-state branching for the reverse reaction experi- 
ments. The resulting rate factors Na < crv > are nor- 
malized to the reference factor Na < crv > re f. determined 
from the reference cross section tr re f., i.e. the MAT data 
with all revisions shown in Table [TT1 The results are shown 



G 



in Fig. [2j For comparison, a theoretical prediction using 
the statistical model is also shown (l7j . 

It has to be pointed out here that a comparison be- 
tween the average cross sections a in SAL and the res- 
onance data in MAT and GRO is possible for the first 
time (see also Fig. [TJ because the total widths V of the 
states in 22 Mg are now available from Table [TJ 
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FIG. 2: (Color online) Ratio between the reaction rate fac- 
tors Na < ctv > from different studies normalized to the 
reference rate factor Na < crv > re f. from the MAT data 
and the modifications in Table [TT1 Long-dashed magenta line: 
5 adopted resonances of the GRO experiment; green short- 
dashed line: original MAT data; black dotted line: MAT data 
with calculated strengths ury for the 5 resonances of GRO 
(MAT-th); horizontal blue dashed and dash-dotted lines: av- 
erage of the SAL data and average multiplied by a factor of 
three to correct for the ground-state branching. (Color codes 
and line styles are identical to Fig. [1]) The conversion of the 
experimental SAL and ANL data to the shown Na < <jv > 
data points is explained in the text. The recommended rate 
factor Na < crv > recommended is located in the narrow over- 
lap of the error bars of the reverse reaction data (SAL x3, 
lightgrey shaded) and the reference rate (grey shaded) at ap- 
prox. 0.55 x Na < ov > rc f. (dark grey shaded). Numerical 
values are listed in Table fVTl A theoretical prediction in the 
statistical model is also shown (brown dash-dotted) [l7j . 



Are the transfer data compatible with the 
direct data? 



Here a strict comparison of experimental data is lim- 
ited to excitation energies E* or resonance energies E 
which have been determined in the MAT transfer and the 
GRO direct experiment. 5 of 8 resonances measured by 
GRO have been uniquely assigned to states which have 
been seen in the transfer experiment of MAT. Two fur- 
ther resonances have not been seen in the MAT transfer 
experiment, but have been detected in other transfer ex- 



periments (see Table V in MAT). The lowest resonance 
in GRO is tentatively assigned to a doublett of states in 
the MAT data. Thus, there is no strict contradiction be- 
tween the experimental transfer data and the direct data 
by GRO. 

For the calculation of reaction rate factors Na < <jv > 
from Eq. JT]) MAT have used the strengths of the 5 
uniquely assigned resonances. The remaining 3 reso- 
nances of GRO are neglected. However, some other states 
have been observed by MAT in the energy range of the 
GRO experiment, and theoretical strengths have been 
used in MAT for these resonances. Consequently, the 
calculated N A < crv > of MAT and GRO agree well 
within the temperature range where the observed reso- 
nances in GRO define the rate factor Na < <jv > whereas 
the Na < <jv > oi GRO is much smaller at lower tem- 
peratures (see Fig. 

A further compatibility test can be made. The reso- 
nance strength uj"f has a general upper limit: uj^/ < loT/2. 
The total widths V v ary between 21 and 144 keV for the 
resonances m GRO hj. These widths impose an ex- 
perimental upper limit on the rate factor Na < crv > 
which is about one order of magnitude above the GRO 
result. The relatively large strengths in GRO are thus 
fully compatible, but at the upper limit of the allowed 
range, especially in the expected case of r a <C T p sa T. 

This finding is further strengthened by a theoretical 
calculation for the only state with a firm spin assign- 
ment. The 2 + state at E — 2130 keV has been considered 
as a member of a Q = 10 higher-nodal rotational band 
in a simple cv-cluster model of 22 Mg = 18 Ne <g> a with 
the semi-magic (N ~ 8) 18 Ne core (similar to [H, (H). 
Even for well-established a-cluster states in the neigh- 
boring nuclei 19 F = 15 N <g> a and 20 Ne = 16 g> a with 
their semi-magic (N = 8) or doubly-magic (N = Z = 8) 
cores it has been found that the calculated width F a 
in this model overestimates the experimental width by 
at least 30 % and often by about a factor of two. Here 
we find for the 2 + resonance at 2130 keV an experimen- 
tal strength of W7 = lO^l 8 '^ keV which corresponds to 
T a = 2.061^;^ keV for T a «C T p . The cluster model pre- 
dicts T a = 1.86 keV. Again, the experimental results of 
GRO are at the upper limit of the allowed range. 

The same procedure has been repeated for the other 
4 resonances from GRO which had entered the original 
calculation of the rate factor in MAT. The results are 
listed in Table ITVl In all cases the a width derived from 
the resonance strengths of the GRO experiment are close 
or even above the theoretical upper limit. This is a very 
unusual finding. However, because only tentative spin as- 
signments are available for these remaining 4 resonances, 
it is not possible to strictly exclude the GRO results from 
the above theoretical considerations. This holds in par- 
ticular for the [4 + ] resonance at 2.287 MeV where the 
GRO width exceeds the theoretical limit by almost a fac- 
tor of 10. In this case the reduction factors from Table 
IIIII suggest J < 2 to allow for the observed resonance 
strength of GRO. 
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TABLE IV: Properties of 5 states in 22 Mg from the GRO data and comparison of the derived from ujT a ~ CJ7 to a 

theoretically estimated maximum width r*' max (discussion see text). All energies E are given in MeV; widths F and resonance 
strengths 017 are given in keV. 

£><JRO £.*t;RO pUHO £;MAT ^Srrorr — pMAT JJ ^ URO pUKO pth.max 

2.17±0.14 10.312 130±80 2.130 10.272 21±3 2* 1031^1 2mf%]£ 1.87 

2.28±0.15 10.422 210±100 2.287 10.429 144±26 [4+] 7.3±?;£ 0.81±J;?? 0.093 

2.52±0.14 10.662 100±50 2.509 10.651 73±19 [3~] 18.2±f"| 2.60± '27 L81 

2.72±0.14 10.862 210±10 2.731 10.873 40±12 [0+] 45.2±^'j 45.2±^;| 97.1 

2.87±0.14 11.012 100±20 2.859 11.001 136±13 [4+] 8.1±| ;g 0.90±g;^ 1.23 



B. Are the transfer data compatible with the 
reverse reaction data? 

A direct comparison between the MAT excitation en- 
ergies and resonance strengths and the SAL average cross 
sections is difficult because there are no common observ- 
ables in the different experimental approaches. Never- 
theless, a comparison can be made in the following way. 
From Eq. (g} and <7 re f.(I?), see Sect. Ill A 11 the average 
cross section a(a,p) can be calculated for each data point 
of SAL and ANL. These calculated a(a,p) should be 
about a factor of m 3 larger because the reverse reac- 
tion data determine only the ground-state contribution 
(a,po). The factor of 3 is taken from the HF calcu- 
lations in SAL (see also discussion in Sect. Ill C"Tj) . The 
results are listed in Tabic [V] 



TABLE V: Average cross sections a(a,p) from a Te f.(E) com- 
pared to the SAL and ANL data. Further discussion see text. 



E cS (a,p) 


exponent 


0"cxp 


Ref. 


Cref. 


Cexp/CVef. 


(MeV) 


for a in mb 










1.194±0.130 


10" 4 


5.5+StS 


SAL 


37.8 




u.i4o_ 067 


1.379±0.129 


10 -3 




SAL 


61.6 


0.062t°;^ 


1.683±0.121 
1.758±0.069 


10~ 2 
10" 2 


n q+0.7 
*' J — 0.8 

3.1 ±0.6 


SAL 
SAL 


33.7 
25.1 


0.068«° 2 5 
0.124 ±0.024 


1.970±0.117 


10- 1 


1 7 +o.7 
^■'-o.e 


SAL 


9.7 


? +0.073 
u - 1 ' d -0.062 


2.568±0.061 


10° 


1.7 ±0.3 


SAL 


20.7 0.082 ±0.015 


1.748±0.077 


10~ 2 


7 c+3.8 
' -°-3.3 


ANL 


22.0 


u.oy:i_ 150 


2.551±0.077 


10° 


2 0+ 10 

z - u -0.9 


ANL 


17.8 


u - 111 -0.053 



The ratios a exp /a rc {. vary between 0.06 and 0.15 for the 
SAL data with a geometric mean of 0.101. The higher 
energy data point of the ANL data is in good agreement 
with the corresponding SAL data point with a ratio of 
0.11 compared to 0.08 from SAL. The lower point of ANL 
is slightly higher by a factor of about 2.5 but has an 
uncertainty of a factor of two. The upper limits of the 
ANL data are compatible with the SAL data points at 
lower energies. Thus, in general the experimental data 
of SAL and ANL are in reasonable agreement. 

As pointed out above (see Sec. Ill C|) . a calculation of 
reaction rate factors Na < crv > from experimental cross 
sections with relatively large uncertainties in the energy 
may have large uncertainties and thus may be misleading. 
Here we estimate the rate factors from the reverse reac- 
tion experiments in the following way. We adopt the en- 



ergy dependence of cr Te {.(E) and use the ratio cr exp /a le {,: 
N A < crv > cxp (T) = N A < av > rof . (T) x !° xp {^ (5) 

cr rc {\E) 

The temperature T for each data point is taken from 
the most effective energy E c g of the Gamow window 
which is given by the well-known relation E c g/kcV = 
122 x (Z^Z^ArodT^) 1 / 3 . This leads to the experimental 
data points shown in Fig. [2] Their average value is al- 
most exactly a factor of 10 lower than the reference rate 
Na < crv > rc f. (horizontal blue dashed line in Fig. [2]). Be- 
cause of the experimental uncertainties of the SAL and 
ANL data, a better determination of the temperature 
dependence of the rate factor Na < ov > from experi- 
mental data for the reverse reaction is not possible. 

A strict comparison of <7 re f. and a eX p requires two 
theoretical considerations. First, calculated resonance 
strengths uj^i have been used in the calculation of 
a le {.(E), and, second, the calculated ground-state 
branching of about one third (as suggested in SAL, see 
also Sect. Ill C 1|) is responsible for an expected factor of 
three discrepancy between <7 rG f. and cr eX p. However, the 
ratio CTexp/tfref. turns out to be about 0.1 (see Table IVf ; 
i.e., it is a further factor of three smaller than expected. 
Both calculations (tvy for the transfer data, the ground 
state branching for the reverse reaction data) are based 
on simple but reasonable arguments, and the uncertain- 
ties should not exceed a factor of two. This factor of two 
for the uncertainties of Na < ov > rc f. and N a < crv > 
from the reverse reaction data are shown as shaded areas 
in Fig. [2] Thus, the real rate can be estimated as follows. 

From the transfer data the reference rate Na < 
av > re f. is derived with an uncertainty of a factor of 
two; i.e., the real rate should be located in the interval 
between 0.5 x Na < crv > Te {. and 2.0 x Na < crv > re f.. 
From the reverse reaction the best estimate for the real 
rate is 0.3 x Na < crv > Te f. (taking into account the 
ground state branching of about 1/3 as discussed above), 
again with an uncertainty of about a factor of two; i.e., 
the real rate should be located in the interval between 
0.15 x Na < crv > re f. and 0.6 x Na < crv > re f.. Combin- 
ing the above intervals, for the real rate only a narrow 
window around (0.5 — 0.6) x Na < crv > re f. remains to 
be compatible with all experimental results. Uncertain- 
ties for this finally recommended rate will be given in 
Sect. [IV] 
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C. Are the direct data compatible with the reverse 
reaction data? 

A strict comparison between the direct GRO data and 
the reverse reaction data is possible because GRO have 
determined the main decay branch of the resonances seen 
in their experiment. In particular, according to GRO, a 
resonance at 2.52 ± 0.14MeV with a width of T = 100 ± 
50 kcV decays mainly to the po channel with a resonance 
strength of = 18.2±f-| keV. MAT have assigned J* = 
[3~], E = 2.509 ± 0.013 MeV, and T = 73 ± 19keV. The 
experimental data points at the highest energies of the 
SAL and ANL data around 2.5 MeV are affected by this 
resonance. Using the resonance energy E from MAT, 
the total width from the MAT experiment (see Table 
H|, and the resonance strength cwy from GRO, we find 
a = 40.2 mb for the SAL data point at 2.568 MeV which 
is a factor of 24 higher than the experimental value of 
1.7 ± 0.3 mb. A similar factor of 23 is found for the ANL 
data point of 2.0+J'gmb where a = 45.3 mb is calculated. 
A much better agreement is found if the huge resonance 
strength of w-f = 18.2 keV from the GRO data is replaced 
by the calculated strength of cwf = 1.12 keV (see Table 
P): a = 2.47 mb for the SAL data point and 2.79 mb for 
the ANL data point. The role of this resonance turns 
out to be minor for the total reaction cross section (see 
Fig. [1]) which is dominated by the long tail of the strong 
0+ resonance at 2.731 MeV. 

This leads to the clear conclusion that there is a strict 
contradiction between the experimental data of GRO in 
the direct experiment and the experimental SAL data us- 
ing the reverse reaction. This conclusion is independent 
of any theoretical calculations. It seems more likely that 
there is an experimental problem in the normalization 
of the GRO data because the unpublished ANL data for 
the reverse reaction are in agreement with the SAL data, 
and thus there must be a problem in the two independent 
SAL and ANL experiments if the GRO data are correct. 



IV. RECOMMENDED REACTION RATE 
FACTOR N A < crv > 

The above discussion leads to the following recom- 
mendations for the reaction rate factor Na < crv > of 
the 18 Ne(a ,p) 21 Na reaction. The most realistic estimate 
from the overlap of the uncertainties in Fig. [2] is located 
around 0.55 x Na < crv > ro f.- Consequently, this rate 
factor is recommended for further use in astrophysical 
calculations. Numerical values are listed in Table IVT1 

Uncertainties for Na < crv > may be estimated as 
follows. A realistic lower limit can be taken from the 
SAL data (multiplied by a factor of three to take into 
account the ground-state branching) which is shown as a 
blue dash-dotted line in Fig. [2j A realistic upper limit is 
the reference rate factor Na < crv > re f.. A strict lower 
limit provide the SAL data (without the correction of the 
ground-state branching, see blue dashed line in Fig. [2]). 



TABLE VI: Recommended rate factor Na < crv > recommended 
in cm 3 s _1 mol _1 and realistic lower and upper limits. Note 

that Na < crv >rec om meiided = 0.55 x Na < crv > re f. . The 

lower limit is given by 0.30 x Na < crv > ro f., and the upper 

limit is given by Na < crv > re f.. 

To, recommended lower upper 

In 4.1xl0 -24 2.3xl0 -24 7.5xl0 -24 

0.2 1.7xl0 -15 9.6xl0~ 16 3.2xl0~ 15 

0.3 3.0xl0" n 1.6xl0~ n 5.4xl0~ n 

0.4 l.lxlO -08 6.1xl0~° 9 2.0xl0~° 8 

0.5 7.2xl0~° 7 4.0xl0" 07 1.3xl0" 06 

0.6 1.7xl0~° 5 9.3xl0~° 6 3.1xl0~° 5 

0.7 2.1xl0" 04 l.lxlO" 04 3.8xl0" 04 

0.8 1.7xl0~° 3 9.3xl0~° 4 3.1xl0~° 3 

0.9 l.OxHT 02 5.7xl0" 03 1.9xl0" 02 

1.0 4.8xl0~° 2 2.7xl0~° 2 8.8xl0~° 2 

1.1 1.8xl0~ 01 1.0xl0~ 01 3.3xl0" 01 

1.2 5.7xl0~ 01 3.2xl0 -01 1.0x10+°° 

1.3 1.5xl0 +0 ° 8.5xl0~ 01 2.8x10+°° 

1.4 3.7x10+°° 2.0x10+°° 6.7x10+°° 

1.5 8.0x10+°° 4.4x10+°° 1.5xl0+ 01 

1.6 1.6xl0+ 01 8.9x10+°° 2.9xl0+ 01 

1.7 3.0xl0+ 01 1.7xl0+ 01 5.5xl0+ 01 

1.8 5.4xl0+ 01 3.0xl0+ 01 9.8xl0+ 01 

1.9 9.1xl0+ 01 5.0xl0+ 01 1.7xl0+° 2 

2.0 1.5xl0+° 2 8.3xl0+ 01 2.7xl0+° 2 

2.1 2.4xl0+° 2 1.3xl0+° 2 4.3xl0+° 2 

2.2 3.7xl0+° 2 2.0xl0+° 2 6.7xl0+ 02 

2.3 5.6xl0+° 2 3.1xl0+° 2 1.0xl0+° 3 

2.4 8.2xl0+° 2 4.5xl0+° 2 1.5xl0+° 3 

2.5 1.2xl0+° 3 6.6xl0+° 2 2.2xl0+° 3 

2.6 1.7xl0+° 3 9.3xl0+° 2 3.1xl0+ 03 

2.7 2.4xl0+° 3 1.3xl0+° 3 4.3xl0+° 3 

2.8 3.3xl0+° 3 1.8xl0+° 3 6.0xl0+° 3 

2.9 4.5xl0+° 3 2.5xl0+° 3 8.1xl0+° 3 
3.0 6.0xl0+° 3 3.3xl0+° 3 l.lxl0+ 04 



A strict upper limit in the astrophyically most relevant 
temperature range is about a factor of three higher than 
the reference rate factor. It can be taken from the GRO 
results; although these results seem to be questionable, 
it has been shown that the GRO strengths are close to 
theoretical upper limits and thus suitable to provide an 
upper limit. 

Finally, this leads to a recommended reaction rate fac- 
tor N A < CTV > recommended = 0.55 X N A < CTV > re f. with 

a realistic uncertainty of a factor of 1.8 and an extreme 
uncertainty of a factor of 5.5; it is interesting to note 
that after all the estimates of uncertainties one ends up 
with an almost Gaussian uncertainty distribution with a 
factor of 1.8 for a realistic (1 sigma) uncertainty and a 
factor of about 5.5 for an extreme uncertainty (3 sigma). 

Of course, the reduction factor of 0.55 between the ref- 
erence calculations and the final recommended rate factor 
also has to be applied to the reference cross section cr r ef. 
shown in Fig. [TJ The absolute value of cr r ef . depends on 
the calculated resonance strengths ljj which are propor- 
tional to the decay widths T a into the a channel (for 
r Q < T p ). Thus, all calculated T a of MAT and in Table 



4.1x 
1.7x 
3.0x 
l.lx 
7.2x 
1.7x 
2.1x 
1.7x 
l.Ox 
4.8 x 
1.8x 
5.7x 
1.5x 
3.7x 
8.0x 
1.6x 
3.0x 
5.4x 
9.1x 
1.5x 
2.4x 
3.7x 
5.6x 
8.2x 
1.2x 
1.7x 
2.4x 
3.3x 
4.5x 
6.0x 



10" 24 

io- 15 
io- 11 

1Q -08 
10 -07 
10 -05 
1Q -04 
1Q -03 

IO" 02 

10 -02 

10 -oi 
10 -oi 

10+°° 
10+°° 
10+°° 
10+ 01 
10+ 01 
10+ 01 
10+ 01 
10+ 02 
10+ 02 
10+ 02 
10+ 02 
10+ 02 
10+ 03 
10+ 03 
10+ 03 
10+ 03 
10+ 03 

10 +03 



2.3xl0" 24 
9.6xl0 -16 
1.6xl0" n 
6.1xHT° 9 
4.0xl0"° 7 
9.3xl0~° 6 
l.lxlO" 04 
9.3xl0~° 4 
5.7xl0"° 3 
2.7xl0~° 2 
1.0xl0~ 01 
3.2xl0~ 01 
8.5xl0" 01 
2.0x10+°° 
4.4x10+°° 
8.9x10+°° 
1.7xl0+ 01 
3.0xl0+ 01 
5.0xl0+ 01 
8.3xl0+ 01 
1.3xl0+° 2 
2.0xl0+° 2 
3.1xl0+° 2 
4.5xl0+° 2 
6.6xl0+° 2 
9.3xl0+° 2 
1.3xl0+° 3 
1.8xl0+° 3 
2.5xl0+° 3 
3.3xl0+° 3 



7.5xl0" 24 
3.2xl0~ 15 
5.4xl0~ n 
2.0xl0" 08 
1.3xl0"° 6 
3.1xl0 -05 
3.8xl0"° 4 
3.1xl0~° 3 
1.9xl0"° 2 
8.8xl0~° 2 
3.3xl0" 01 
1.0x10+°° 
2.8x10+°° 
6.7x10+°° 
1.5xl0+ 01 
2.9xl0+ 01 
5.5xl0+ 01 
9.8xl0+ 01 
1.7xl0+ 02 
2.7xl0+ 02 
4.3xl0+ 02 
6.7xl0+ 02 
1.0xl0+° 3 
1.5xl0+° 3 
2.2xl0+° 3 
3.1xl0+° 3 
4.3xl0+° 3 
6.0xl0+ 03 
8.1xl0+° 3 
l.lxlO+° 4 
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[II] should be reduced by the same factor of 0.55. 

The recommended rate factor Na < crv > recommended 
is slightly lower than the MAT rate factor around T9 w 1, 
but significantly smaller at higher temperatures around 
T 9 « 2. The minor difference between the MAT rate 
factor and Na < crv > re commended a t lower temperatures 
is due to a compensation of the enhancement from new 
spin assignments from the CHA transfer data (leading to 
smaller spins J and thus increased resonance strengths) 
and the derived reduction factor of 0.55 from the com- 
parison with the reverse reaction data. The significant 
decrease at higher temperatures is mainly a consequence 
of the replacement of the huge resonance strengths from 
the GRO experiment by smaller calculated resonance 
strengths. 

In a comparison of the new recommended rate factor 
Na < crv > r ecommended with the SAL result (sec Ta- 
ble II in J{|) it has to be kept in mind that SAL pro- 
vide the rate factor for the ground-state contribution 
18 Ne(a,po) 21 Na g . s . which is derived from their experi- 
mental data; the given upper and lower limits are calcu- 
lated from their experimental uncertainties, but do not 
include the additional contributions of the {a,Pi^o) chan- 
nels and the corresponding uncertainties. Thus, it is not 
surprising that the new recommended rate factor exceeds 
the upper limit of the SAL rate factor. The new recom- 
mended rate factor is about a factor of 5 higher than the 
SAL result in the astrophysically most relevant temper- 
ature range of T 9 = 1 — 2. 

The theoretical prediction of the 18 Ne(a,p) 21 Na cross 
section in [17J is based on the statistical model and a 
global parameter set. The calculation has been done 
before the experimental results of GRO, SAL, MAT, 
ANL, CHA, and MAT were available. By definition, 
such a statistical model calculation cannot reproduce 
details of the <j(E) curve shown in Fig. [TJ Neverthe- 
less, the predicted rate factor Na < crv > is in reason- 
able agreement with the recommended rate factor (see 
Fig. [2]) and remains within the realistic uncertainty esti- 
mate of Na < crv > recommended m the temperature range 
Tg = 1 — 3. However, the temperature dependence of the 
theoretical rate factor is slightly steeper compared to the 
recommended rate factor, and thus at low temperatures 
below T 9 «l the theoretical rate factor is located below 
the recommended rate factor between the realistic and 
the the extreme lower limit of the recommendation. 

The recommended rate factor Na < crv > recommended 
is fitted by the usual expression, see e.g. [T3|, Eq. (16): 



N A < crv > 
cm 3 s~ 1 mol _1 



exp (a + aiT 9 1 + a 2 T g 1/3 + a 3 Tg 



-,1/3 



-a 4 T 9 + a 5 T^ 3 + a 6 In Tg) 



(6) 



The a,i parameters are listed in Table I VIII The deviation 
of the fitted rate factor is always below 10 % over the full 
temperature range 0.25 < Tg < 3 and typically below 
5 % in the most relevant range 1 < Tg < 2. 



V. SUMMARY AND CONCLUSIONS 

The present knowledge of the reaction rate factor 
Na < crv > of the 18 Ne(a,p) 21 Na reaction has been sum- 
marized. For this purpose experimental results from dif- 
ferent experimental techniques are combined. Transfer 
reactions provide the best determination of the excita- 
tion energy E* and spin and parity J* of states in 22 Mg 
which appear as resonances in the 18 Ne(a,p) 21 Na reac- 
tion; however, transfer reactions cannot provide the re- 
quired resonance strengths wy. These strengths have to 
be taken from theory or from direct experiments which 
are however extremely difficult and require the combina- 
tion of a radioactive 18 Ne beam and a helium gas target. 
Complementary information has been derived from the 
experimental study of the reverse 21 Na(p,a) 18 Ne reaction 
using a radioactive 21 Na beam and a solid CH2 target. 

A basic prerequisite for the comparison of results from 
various experimental techniques is the availability of to- 
tal widths T for the resonances under study. The total 
widths r were determined from a rcanalysis of the peak 
widths in the MAT experiment. 

A compatibility test between the results from various 
experimental techniques shows that there is no contra- 
diction between the various experimental data except the 
disagreement between the direct GRO data and the re- 
verse reaction data from SAL and ANL. This leads to the 
conclusion that the most likely explanation is a problem 
in the normalization of the GRO data. Consequently, 
resonance strengths from GRO have been replaced by 
theoretical resonance strengths in the calculation of the 
rate factor Na < crv >. 

The calculation of Na < crv > for the 18 Ne(a,p) 21 Na 
reaction from transfer data requires theoretical resonance 
strengths, and the calculation of Na < crv > from the re- 
verse 21 Na(p,a) 18 Ne reaction data requires a theoretical 
estimate of the (a,po) ground-state branching. Both cal- 
culations are based on simple but reasonable arguments, 
and the corresponding uncertainties should not exceed 
a factor of two. This leads to a relatively narrow over- 
lap region between the higher Na < crv > calculated 
from transfer and the lower Na < crv > calculated from 
the reverse reaction data. This narrow overlap region is 
considered as the new recommended reaction rate factor 
Na < crv >recommcnded- The uncertainty of the recom- 
mended rate factor is about a factor of 1.8 (sa la uncer- 
tainty). For Tg = 1 — 3 a theoretical prediction [ijj lies 
within this error band, but the theoretical temperature 
dependence of the rate factor Na < crv > is somewhat 
steeper than the new recommendation. 

The new recommended rate factor is slightly lower 
than the MAT rate factor at low temperatures and sig- 
nificantly smaller at higher temperatures, and the new 
rate factor exceeds the SAL result by about a factor of 5. 
The strong conclusion of SAL (based on their lower limit 
for the rate factor Na < crv >) that "the breakout from 
the HCNO cycle via the 18 Ne(a,p) 21 Na reaction is de- 
layed and occurs at higher temperatures than previously 
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TABLE VII: Fit parameters a,i of the reaction rate factor Na < crv > from Eq. ©. 



do 


ai 


a.2 az 


a 4 a 5 


a 6 


-21.0595 


-0.3301 


-58.1167 89.3359 


-14.8713 1.9862 


-24.1080 



predicted" cannot be supported. Instead, because of the 
only minor deviations of Na < crv > re commended from the 
MAT result at low temperatures around Tg = 1, the ear- 
lier conclusions of MAT should remain valid in general. 
Further astrophysical network calculations with the new 
recommended rate factor Na < crv > rC commcndcd are re- 
quired to study the relevance of the modified temperature 
dependence of the rate factor in detail. 
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